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Abstract 

It is shown that quantum devices based only on oscillators cannot serve 
as the universal quantum computer, despite of entanglement in such devices, 
'. which we roughly estimate for the ideal case and for the harmful entanglement 

with photonic modes. We show that quasi-particles are the native shell for the 
entanglement already for ground state, in contast to the free electromagnetic 
field where vacuum state does not produce entanglement at all. 

1 Introduction 

Entanglement of a wave function ^(ri, r 2 , . . . , r n ) of n particles means that it cannot 
be represented as the product ^i(r 1 )^ 2 ( r 2) ■ ■ ■ ^n( r n) of one particle wave functions. 
The analogous definition takes place for the entanglement between two subsets of 
{1,2, ... ,n} that is the so-called bi-partite entanglement. This fundamental no- 
tion plays the key role in quantum information processing, especially in quantum 
computing, where entanglement is the necessary property of quantum states, which 
arise in the course of quantum computations, if it pretends to be faster than its 
classical counterpart. Entanglement is the physical phenomenon, which results in 
quantum non locality, that is the violation of Bell inequality, detected in numerous 
experiments. 

The important role, which entanglement plays in information processing explains 
that the most elaborated numerical measures of this property used the discrete form 
of quantum state representation, e.g. its qubit form 

N-l 
3=0 
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where N = 2 n , n is the number of qubits. In this representation the product of wave 
functions turns to be tensor product, an we can pass to the continuous representation 
and vice versa with some lack in the accuracy that follows from the approximate 
equation 

= 

j 

where 5j is proportional to the characteristic function of j-th cube in the configu- 
ration space for the continuous wave function, which plays the role of Because 
of exponentially huge dimentionality of Hilbert space (values of J) the last equa- 
tion cannot be accurate even for relatively small n, and just this non accuracy 
dramatically complicates the application of discrete entanglement measures to the 
continuous case. 

The simple and robust measure of entanglement for many qubits (and qu-dits 
- d- levels quantum elements) has been proposed in the paper [5]. It is the min- 
imal value Ch(^>) of Shannon entropy of quantum state amplitude distribution: 
~~ IAj'I^GAj'I) taken after all possible one-qubit unitary operations on the given 
quantum state Entanglement measure Ch possesses all properties of entropy, 
including the different form of additivity, and has extremal value on W- and GHZ- 
states; there are efficient algorithms for its computation and it can be generalized 
on fermionic quantum states. But this measure hardly can be applied to states of 
continuous objects, like the field, because of its substantially discrete character. 

Instead of it we use the simple idea (look at [1]) how to check entanglement for a 
function with two argument by the sequential differentiation. The function f(x,y) 
is non entangled if and only if 

d 2 ln f n 

= 0. 

dxdy 

This will be the indicator of bi-partite entanglement if we sum up this values over all 
x G X, y EY for the division X, Y of our configuration space. Following Born rule, 
we obtain the bi-partite measure of entanglement for the continuous wave function 
*&(ui, U2, ■ ■ • , u n ) in the form 

Eij = / 1* 

K 

The measure of bi-partite entanglement given by (1) has the transarent sense 
in the case of kernel- type wave functions of systems with quadratic Lagranjians, 
which have the form ^ = F(t)exp(^) where S is the action along the classical path. 
For example, for the interaction between oscillating charged particles and the field, 
represented by the amplitude a the action has the form S = S part i C i es + Sf ie u — \ ja- 
For the small t the value of our measure will be independent of the coordinates: 
Eparticies, field = j^- For the general wave functions it is not true, because the 
summing of such values on the different paths changes entanglement. We will use 



d 2 ln g 
dmduj 



du\du2 ■ ■ ■ du n 



(1) 



2 



this measure for the entanglement of the two kind of media: interacting harmonic 
oscillators and electromagnetic field. 

However, in some cases it is more convenient to use the canonic expression of 
the bi-partite entanglement measure, given by the formula 

Mij = -tr(pi ln(pi)) = -tr(pj ln(pj)). (2) 



2 Interacting oscillators 

A system of interacting oscillators form the basis for the most elaborated technology 
of quantum computer: excited levels of ions in Paul trap. There quantum gates on 
excited (|1)) and ground state (|0)) levels of ions are organized with the ancillary 
qubits, which role plays phonons of the mechanical oscillations of ions in trap (see 
[10], [3]). For the entangling of ion levels with the oscillation mode classical laser 
field is used. This model: real charged oscillators plus field form the abstract model 
of such a quantum computer. The field can be represented as the ensemble of 
(imaginary) oscillators as well; we factually will take up entanglement in oscillator 
systems. We start with the one dimensional system of real charged particles of the 
same mass, each of which has its equilibrium position, and let Uj denote the shift 
from the equilibrium for j-th particle. If k is Gook coefficient for the force between 
the neigbor particles, Hamiltonian of the whole ensemble has the form 



N 



'2m 



n=l 



Let the total number N of oscillators be so large that we can ignore boundary 
conditions (alternatively, we can enclose the string of oscillators to the ring without 
bounds). When passing to a media (field) we launch N to infinity and use integrals 
instead of sums. 

Even for not large N it is impossible to operate with this Hamiltonian directly 
because of the summand of interaction u n u n+ i, which cause entanglement. We must 
use the canonic transformation to get rid of interaction: 



7; — — V II P - i( i nd 

q (A) 

TJ — _J_ P iqnd \ 1 



n 



where d = 2ir/N. 

This is the canonic transformation, which elimitanes the interaction; in new 
coordinates U we have the system of non interacting oscillators. 

Using the definition of impulse p n = j^- and rules of differentiation, we obtain 
formulas for the transformation of impulses: 
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_ _j_ p p iqnd 

<? (K) 
p — _J_ n p -iqnd \ U J 



where P n = -7J7- 

n • dU„ 



Canonic transformation is linear and transforms any small cube of the division of 
the configuration space to the cube due the orthogonality of (4) . We thus can treat 
the canonic transform as the permutation of basic vectors in Hilbert space of quan- 
tum states, like CNOT two-qubit operation: \x,y) — > \x,x + y (mod 2)). Canonic 
trasformation is one-to-one corresponding between points in the configuration space 
for our system of oscillators (see Figure 1). This transformation is entangling and 
we will estimate the measure of entanglement it produces. It is important, that this 
transformation eliminate the initial entanglement, which means that it incapsulate 
this entanglement into quasi-particles, called phonons. Spectrum of phonons is dis- 
crete (see (6)); we will see that the discrete form of spectrum is typical for ensembles 
of oscillator type with entanglement. 

Turn the coordinate system for q so that this parameter, replacing n, takes value 
from the symmetric segment. Instead of q + q' = N we then write q + q' = 0. The 
relation of inequality is induced form the old set 1,2, ... ,N so that pairs where 
q > —q are approximately the half. 

Rewriting Hamiltonian in the new coordinates, we have: 

H= E 2^v(E P q P q >e^'> d ) + § E U q U q , 

n=l q,q' q 

-%^(U q U q ,e- iqnd e- iq '( n+ V d ) = 
q,q' 

N 

= E PqP-q " f E E e^^e"*' + f E U q U q , = 

q q,q' n=i q 

-J-VPP -KSTTJTJ P +iqd 1 K tj tj 
2mN £-1 1 1 N ^-t 1 1 ~ N ^ 1 1 

= 2^vE PqP-q + f ^ U q U- q (l - COs\qd)). 

q q>-q 



Here K = mu 2 /2, and the standard formula for summing of the geometric pro- 
gression of exponents was used, which gives zero for q 7^ q' , and also we put pairs 
q, —q into order so that explicitly written is only a half: for which q > q'- that gives 
coefficient 2 in the last summand. 

We pass again to the new variables, real numbers: 

77 — Y -L- lV Y — U 1 + U -1 V — Uq-U-q . 

V q — SL q ~rll q , s^ q — 2 , l q — 2i , 

x q = -^Y, u nCos(qnd), Y q = -j- u n sin(qnd) , 

n n 
d . _d_l 1 d 1 

dU q dX„2 dY„ 2V 
d _ a 1 _ 8 1 
dU- q dX q 2 dY q 2i ' 

d 2 __ 1 ( d 2 , d 2 
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Figure 1: Quasi-particles 
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A last obtain 

ff =-ii:(44) + ^ {x < + 17,(1 ^ cos(9d)) ' 

We see that in the new coordinates our system represents the set of independent 
harmonic oscillators of the mass m = 2m, with new coefficient K = 2K(1 — cos(qd)) 
and frequences 

j2K 

0J q =\ (l-cos(qd) (6) 

V m 

It means that the ground state of these new oscillators, which are quasi-particles, 
is the product of ground states of the form 

*„(*) = -j= [-^) exp (-_ j H^muj/h) 

(x- amplitude, n — 0) of each of them separately (the same is true for excited states: 
they are excited independently of each other). Since Hermitian polynamials have 
the form H — 1, Hi — x, H 2 — x 2 — 1, the ground state (with H ) can be written as 

^ne)" ! -(-=«). 



Transforming this into old coordinates we have 



The expression inside square brackets is U q U_ q , and if all u q be equal it would not be 
entangled state, because U q U_ q = v? n . However, accordingly to (6) all u q are 

q n 

different and this reasoning is not right. The ground state of the system of dependent 
harmonic oscillators on the level of their initial coordinates turns entangled. 

We estimate the measure of entanglement by the formula (1). It gives the integral 



Eij = J \^o\ 2 So\duid,U2 ■ ■ ■ du n 



where S = ^ —^^cos((i — j)qd). Since So does not depend on u, we found that 



9 



Eij = S . The ground state is entangled because S > 0. For the first excited state 
with one phonon in s-th cosine mode: 



*i» = II {^h) 1 ' 2 ex P ( ~1Kn ( {Y. u nCOs{qnd)f + u n sin(qnd)) 

q \ \ n n 

-^{T. u nCos(snd)) 



Now the application of the definition (1) after simple transformations gives 
Ei,j(®is) — J \ \^is\ 2 So — \^o\ 2 cos(sid)cos(sjd)\duidu2 ■ ■ ■ du n > 0. 

Here we cannot complete without numerical computation; the best way is Monte 
Carlo method for many dimensional integrals. 

For the case of two excitations we have two possibilities: a) excitations of the 
different modes, or b) excitations of the same mode. For these cases we obtain 
correspondingly: 

E i j(^ 2s p) = j\\^2 S p\ 2 S - -^\^i p \ 2 cos(sid)cos(sjd) - i s \ 2 cos(pid)cos(pjd)\ 
du\du2 ■ ■ ■ du n: 

Ei,j(^2 S ) = J ||^2s| 2 5'o + \^o\ 2 2u i cos(sid)2u j cos(sjd)({J2 cos ( snd )) 2 - 3 )l 

n 

du\du2 ■ ■ ■ du n , 

(8) 

Numerical estimations show that the measure of entanglement becomes more 
complicated with the growth of excitation. 

3 Entanglement in the combined systems 

With the purpose to find the adequate description of quantum state dynamics we 
must go beyond the limits of simple quantum mechanics, and include the electro- 
magnetic field into consideration. 

Electromagnetic field is determined by its potential 

(<M) 

here <f> is the scalar field, A is vector field. Electric and magnetic fields are expressed 
as follows: 

dA 

E = -V0- — , B = V x A. 
at 

If for a given function ip we change vector and scalar field as 

A' = A + V^ 0' = 0-^ 

it will have no physical sense; this transformation is thus called calibration and we 
are free to choose it at our convenience. We further fix Lorenz calibration in order 
to write Maxwell equations briefly: 

UA» = j M , = 0, V M A M = 0. 

where A^ = (0, A x , A y , A z ) and □ = V M V M = j£ - V 2 . 
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Quantization of the system "field + particles" presumes that a quantum state 
of this system has the form ty(A^(R, t), q(t)) where is the field, detremined in 
each point of space-time, q(t) are trajectories of all charged particles. Unitary law 
of evolution for \I/ in non relativistic case (when t we consider as the real time for 
all particles) can be represented in the form of Feynman path integrals: 

*(A /( „, ,>,) = / K(Af in , A m , fc)*(A ( „, - qin )DA lin D^ 

where the kernel K results from the integration along all trajectories 7 : A(t), q{t), t G 
[tin,tfin], leading from the initial classical state A in ,q in to the final classical state 
A fi n , q~f in of our system: 

K = J ex P {~S[ 1 ])D 1 
where the action of "field + particles" has the form 

S[ 1 ] = J2y J \4i\ 2 dt+J jlp(R,t)0(R,t)~j(R,t)A(R,t)]d 3 Rdt+-^- J (E 2 -B 2 )d 3 Rdt 
For electrodynamics the canonic transformation has the form: 



A(R,t) =V^cfa k (ty kR ^, ak{t ) = -^-jA(R,t)e~* kR dR 



(2tt) 3 ' v^c 
ikR_£k _ r 11-D +\ a -ikR 



<p(r, t) =j cf> k (ty kR -^, Mt) = J <t>(R, t) 



e 



j(R,t) = fj k (t)e ikR ^, j k (t) = J j(R,t)e^ kR (9) 

p(R,t) =f Pk (t)e ikR ^, p k (t) = fp(R,t)e-* kR , 

Maxwell equations then can be rewritten in the form d\ k + k 2 c 2 a\ k = y/^njik 
and the same for the second polarization, which means that the field is represented 
as the pairs of complex oscillators (see ([7]). 

What is the entanglement of the electromagnetic field? At first we consider the 
pure field without particles. Its psi-function we can represent as Ylk^ni( a k) where 
\& n is n-th excited level of harmonic oscillator. If we pass to the real amplitudes, k 
and —k values of impulse will be grouped and for this "joint" value of k we have 
four independent real modes of oscillators, whereas is we distinguish k and —k we 
will have two independent modes, which we can distribute between two directions 
of polarisation: 1 and 2. Entanglement of the field we can estimate by the measure 
Eij if we applied it to R- representation of the field. For example, the entanglement 
of the ground state (without photons) is Eij(^/ ) where 

^ = Y[ul /2 exp{-^^[{J2arCos{kRd)) 2 + (J^ a r sin(kRd)) 2 }. 

k Ft Ft 

The direct computation gives 

s «(*o) = E w cosN(i ^ _ Ri)] = °' 
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because the sum is the real part of integral f exp(ikR&)dk, for which integration 
over sphere S 2 : \k\ — k gives zero for any k for Ra = Rj — Ri ^ 0. 

On the other side, the analogous estimation for the field with only one photon 
(the single excited mode) shows non zero entanglement; the field with the single 
photon in each mode shows the entanglement in coordinate space between only 
pairs of the form R, —R. 

The entanglement in R- representation of the wave function has the direct con- 
nection with the entanglement of charged particles, interacting with the field. This 
entanglement shows the potential entanglement between the different charged par- 
ticles, which can absorb or emit the photon. For example, let we be given three 
particles, each of which can absorb a photon and pass from the ground (|0)) to the 
excited state (|1)). We then obtain the state of the form 

a(|100) + |010) + |100>) + 01000). (10) 

Of course, the degree of entanglement between particles will not be equal to Eij, the 
last is only the indicator ("finger print") of the "real" entanglement of the charged 
particles. 



4 Quasi-particle as the shield of entanglement 

Entanglement is the fundamental phenomenon, which was demonstrated in numer- 
ous experiments. Its deep nature follows from non-locality of quantum states, which 
entanglement is sufficienty large. Its first visible application concerns quantum com- 
puters, namely the realization of Grover algorithm ([8]). For this aim we have to 
build quantum states of the form 

i*> = «5>>+« ( n ) 

3 ¥=30 

where a 2 (N — 1) + /3 2 = 1. The demonstration of such states would already mean 
that we can realize Grover search algorithm, which is such evolution of \^/) that 
(3 — >oo. 

Can we realize this algorithm by means of only charged particles of oscillator 
type? 

Lagranjian of the oscillator with the external force / has the form 

L = -u 2 - —u 2 + f(t)u 
and its kernel with the interaction (see [7]) is: 

K ™= j*£^r e < s °^ < 12 » 
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where the classical action of the oscillator is 



SZssic = liftr [cos(uT(u f + Ui ) - 2u fUl 



/ /(*)«» "(* - + S / - t)dt 

(13) 

/ / f\ t )f\ s ) sin ^itf ~ t)sinu(s - ti)ds dt], 
u u 
T = t f -U. 



If the oscillator is the photon mode and external force comes from the charged 
particle K osc = K A turns to the photon kernel for the fixed trajectory of particles 
and external force / = j k 1(2) of photon oscillator - to the current of charge. All 
currents we must sum up ocer all participating particles. 

We now consider the model with the conformation of chain of atoms. For ex- 
ample, it can be a protein of nuclear acid, for which the conformation - the spatial 
location of the chain - is conserved by hydroden bonds between segments of the 
chain. Let the conformation be given by the three dimension vector function of the 
real parameter r of the form q = q(r). We treat that oscillators (electrons) are lo- 
cated along this chain so that their small shifts do not depend on the conformation 
of the whole chain in the space. Electrons are in quantum states as well, but we 
account only the oscillations of their density. The parameter r thus plays the role 
of n - number of oscillator from the model example (see the previous section). We 
denote by u T the shift of the oscillator corresponding to the value r. The density of 
charge (electronic cloud) we denote by p(q). It velocity is u T . 

If the charge e was point wise and was located in the point q, its density p(R) 
would be delta-function p = e5(R — q), and the current created by its movement 
would be equal eq5(R—q). Since Fourier transform of delta-function equals e~ tkq , we 
conclude that in the case of one point wise change, which performs small oscillations 
arond the point q 

U k = e(q)u(q)e- ik «, 3k = e(q)U k = e(q)u(q)e- ikq (14) 
Canonic transform for the system of oscillators looks as follows: 

n 

and the velocity 

n 

From (14) we obtain for the oscillations of the continuous distribution of electronic 
density 

U k = J p(q(T))u T e- ikq( ^dT 

Jk = J p(q(r))u T e- ik ^dr (15) 
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where the dot as usually, means the differentiation on the time. The coordinate q 
depends on the internal coordinate on the chain r, but for the fixed conformation 
of the chain in space q doe not depend on the time t. If we treat the deformation of 
the chain, the dependence q T (t) from the time arises. However, its dynamics is much 
more slow, than the dependence u T (t). The dependence u from the time determins 
phonons in the chain, and at last, the emission of photons that is much more rapid 
process than the change of conformation q T , hence at the work with quasi-particles 
we ignore the dynamics of conformation. 

We now find the Lagranjian of the system: " 3-dimensional particle + field", 
following the analogy with the already considered examples. The current j k = 
Fk/rriphotV^K, k — \k\ — kphot, where m p hot is the Active "mass of the photon". 
From the other side j k = ell k where is the coordinate in k- basis of the sfift of 
particle-oscillator of the mass fh = m part and the coefficient of elasticity K = K part 
respectively to its equilibrium position. 

The general Lagranjian L k = L± k or L 2 k- in dependency of the polarisation of 
the field, which we omit, is The 

Lk= ^ _ + ^ eak u kTriphot 

mpart Jj2 m partU k part j j2 ^ ' 

We note that m p hot gives only the multiplier in the kernel and thus can be ignored. 

By the analogy of one particle Lagranjian (16) we can compose Lagranjian for the 
chain with the field and write the action and the kernel analogously to formulas (12), 
taking into account already found action of oscillators of the field at the movement 
of them along the classical trajectory and the fixed movement of particles: (13). 
Here for the currents and shifts of electronic density the expression (15) should be 
used. 

The result looks as follows. The kernel of the system "chain + field" has the 
form: 



K k gen {t)= I exp(if^(jl-ul par M)K k phot (t)DU k , 

^part t% 

K phot(t) = yj 2n2Tin(kct) eX P( \ S phot i U k ( r )} ) > 

S$Z[Uk(r)] = fX^[co8(kct)(al fm + a\ J 

t 

-a k fi n a k in + 2 rT ph Tkc I m phot V47rj k (b)sin(kc(b - t in ))db ( 17 ) 

tin 

t 

+2 ^M I m photV ^ Jk (b)sin(kc(b - t in ))db 

tin 

t b 

- m 2 h 2 otk 2 c 2 J J jk(b)jk(s)ml hot 4:7rsin(kc(t - b))sin{kc(s - t )ds db] 

tin t%Ti 

We see that there is the complete parallelism of the computation on the pairs of 
the form k, —k, which does not take place in the general case. Peculiarity of our case 
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is that particles are included into threads that makes possible the complete passage 
into k- representation for the field and for particles as well. Quantum dynamics of the 
peculiar case we have just considered can be thus effectively simulated on classical 
computers. On the other hand, the ability to produce of the states (11) means the 
ability to fulfil Grover search algorithm, which cannot be effectively simulated on 
classical computers (see [2], [9]). It proves, that to produce the states of the form 
(11) via quantum processors of the type "chains + field" is impossible. 

Entanglement of oscillator type is shown in experiments with ions in trap (see 
[10]); this is the states of W-type: |100) + |010) + |001) and GHZ-type |000) + |111>. 
The oscillator character of these types of entangled states follows from the possibility 
to obtain them by wave-like process. For example, the single photon can cause 
excitation of only one of ions (W-type); for the example from above (10) it is possible 
to make (3 negligible by using resonators. These types is entanglement, protected 
into quasi-particles of a system of oscillators. We have shown that purely oscillator 
type quantum processor is not sufficient to build the robust quantum computer. 

It is important to obtain the states of the form (11), because they do not belong 
to the oscillator type. The experimental demonstration of such states would be 
the decisive argument for the practical validity of the direct approach to quantum 
computer - via quantum gates (see [6]). 

5 Decoherence and energy gap 

The source of decoherence is the entanglement between the considered system and 
its environment caused by photon emission. In the abstract quantum computing 
the formal way to suppress decoherence is quantum error correcting codes ([4] and 
many others) . For this it is necessary to have the robust quantum devices with few 
tens qubits that requires resistence against decoherence on the physical level. We 
now estimate the degree of decoherence via entanglement with the emitted photon 
of the simgle atom and show that the dependence on its frequency is approximately 
linear. It shows that we can control decoherence by avioding of passages with big 
energy gaps even if real phonots are emitted by a quantum device. 

Let \&j and be the initial and final states of the atom correspondingly, ujif = 
(Ei — Ef)/h be the frequency of the passage. Let the polarization of the emitted 
photon be denoted by e , and the only restiction on its impulse direction is that is 
orthogonal to e. The final state of the atom and electromagnetic field will be 



Here the phase of does not influence to the entanglement if we use (2). For the 
small frequences larger than the atom size) we can use the dipole approxima- 
tion: pk = |Aj/| 2 ~ Wj/(e/ij/) 2 , where faf is the dipole electric momentum of the 
atom ([7]). The entanglement measure between the atom and the electromagetic 





k: \kc\=uJif 
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field is thus the integral on the two dimension sphere 

M if = - J 4:7r 2 uj if (efi) 2 dSe 
s 2 

where S s is the element of square on the surface of S 2 corresponding to e. We 
straightforwardly obtain that 

M) if = -An 2 cu if [A ln(4n 2 tu if ) + B] 

where A = J s2 (eJx) 2 dSe, B = f s2 2(eft,) 2 ln(\efi\)dSe. The second root of the quadratic 
polynamial has no physical sense because it does not belong to the area of validity 
of the dipole approximation. It shows that the dependence of entanglement (and 
the decoherence resulted from it) is proportional to the energy gap corresponding 
to the considered passage. 

6 Conclusion 

We have evaluated the bi-partite entanglement in a system of interacting harmonic 
oscillators, for which all eigenstates beginning from the ground state are entangled, 
and have shown that the entanglement of spatial components of electromagnetic 
field can serve as a "finger print" of the entanglement between real particles, which 
emit the field. Entanglement can be estimated for systems of particles and elec- 
tromagnetic field where it turns linear from the energy gap. We have proved that 
the system of oscillators with the field cannot play the role of quantum computer, 
because its evolution can be effectively simulated on a classical computer. The en- 
tangled states obtained in experiments have the type of oscillator states, whereas 
Grover search algorithm requires states of the other type. Grover states cannot arise 
in a system based only on harmonic oscillators, and the experimental demonstration 
of them represents the significant step on the way to quantum computer. 
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